Discerning Incompressible and Compressible Phases of Cold Atoms in Optical Lattices 
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Experiments with cold atoms trapped in optical lattices offer the potential to realize a variety 
of novel phases but suffer from severe spatial inhomogeneity that can obscure signatures of new 
phases of matter and phase boundaries. We use a high temperature series expansion to show 
that compressibility in the core of a trapped Fermi-Hubbard system is related to measurements 
of changes in double occupancy. This core compressibility filters out edge effects, offering a direct 
probe of compressibility independent of inhomogeneity. A comparison with experiments is made. 

PACS numbers: 03.75.Ss, 71.10.Fd 



The search for stable quantum many-body phases 
forms the basis of quantum condensed matter, quan- 
tum chemistry and elementary particle physics. Stable 
states often arise as a consequence of energy gaps that 
set an energy scale for resilience. Examples of gapped 
condensed matter phases include superconductors and 
Mott insulators. In these phases gaps arise in the pres- 
ence of (or as a result of) interactions to form many-body 
states that resist small perturbations from the environ- 
ment. Gaps play a role in spectacular quantum effects 
in the solid state including zero resistance in supercon- 
ductors or a dramatic rise in resistance in Mott insula- 
tors. The high T c compounds arc believed to harbor both 
states in a possible union [l[ that leads to an enhance- 
ment of the superconducting energy gap and therefore T c 
in a doped correlated Mott insulator. Electronic mod- 
els of the high T c compounds (repulsive two-dimensional 
Fermi-Hubbard models) are believed to have rich phase 
diagrams and may even exhibit this unique type of su- 
perconductivity. Yet the low temperature solution of the 
two-dimensional Hubbard model remains elusive. Recent 
work 0, [H seeks to emulate the Hubbard models in an 
ideal setting, with cold atoms in optical lattices, in a 
search for a superconducting state and other phases. 

Optical trapping and cooling of cold atomic gasses 
allows the preparation of nearly ideal manifestations 
of Bose and Fermi-Hubbard models in the laboratory 
[1, H, H 0| ■ A broad variety of many-body phases have 
been predicted Q. However, few techniques exist for 
experimentally observing new phases and their proper- 
ties. Techniques currently in use include time of flight 
imaging of the momentum distribution (predominately in 
bosonic systems) and noise correlations [9( , optical molec- 
ular spectroscopy of pai r correlation functions and 



Bragg spectroscopy [Tl|. Proposed techniques include 
using edge currents in trapped rotating lattices [12[ or 
Fourier sampling of time-of-flight images to reveal new 
correlation functions 



most experiments, present a major difficulty in realizing 
and observing bulk phases of the Hubbard model with 
optical lattice experiments: e.g., a large portion of the 
sample becomes compressible near the edges even when 
a Mott insulator has formed in the center. Extracting the 
physics of a homogeneous bulk system by separating the 
approximately homogeneous behavior in the flat center 
of the trap from the surrounding inhomogeneous system 
is an intrinsically difficult challen ge in these experiments 
and for comparison to theory 14, HI, 1(| 17, 18, 1^, 2cj |. 
In this Letter we use a high temperature series expan- 



21 , 22|,[23| to directly relate recent experiments [f| to 
(a measure of) "core compressibility" , the compressibil- 
ity of the atoms near the center of the sample, excluding 
the edges. We argue that this measure of the core com- 
pressibility can be used to detect stable, incompressible 
phases and offers a valuable tool for mapping out phase 
diagrams. 

We begin our analysis with the Fermi-Hubbard model 
of cold fermionic atoms in optical lattices for equal 
pseudo-spin populations where the spin index refers to 
different hyperfine levels: 
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However, harmonic trapping potentials, inherent in 



Here the site dependent chemical potential, /Xj = n — 
7R 2 varies with the discrete spatial coordinates, = 
(ixi iy,iz), as a result of the parabolically confining laser 
beam waist and magnetic trapping potentials. The s- 
wave interaction yields an on-site repulsive interaction, 
of strength U , which is tunable via a Feshbach resonance. 
The hopping between nearest neighbor sites, t, changes 
with the lattice depth and is therefore tunable with the 
laser intensity. The number operator ni — nf + n\ mea- 
sures the number of fermions at a site i. 

For our quantitative calculations we use parameters 
relevant to recent experiments. Ref. [f| places two hyper- 
fine species of 40 K in a simple cubic optical lattice with 
low hopping t = 0.05 — 0.2 kHz, a tunable on-site term, 
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U = 0-9kHz and 7 = 0.003-0.005 kHz. Here and in the 
following we work in kHz by setting h = 1 and = 1. 
We also set 7 = 0.00384 kHz, unless otherwise noted, 
and choose all other parameters to make contact with 
Ref. @. Suitably chosen chemical potentials, ji ~ 1 — 7 
kHz, yield total particle numbers, N ~ 10 4 — 10 6 . Esti- 
mates in Ref. [(§ find that the largest unknown, lattice 
temperature, can be kept below U with values as low as 
T w 0.8 kHz and possibly lower. 

To theoretically analyze observable signatures of in- 
comprcssibility in a trapped Hubbard model we use a 
high temperature series expansion of the grand partition 
function Z = Trexp(— (3H) about the atomic, t = 0, 
limit. Such high temperature series expansions have 
a long history [23| and yield exact results for thermo- 
dynamic quantities of the Hubbard model in the limit 
/3t <C 1. Note that all experiments done with fermions 
in optical lattices currently lie in this high temperature 
regime when parameters are tuned to U 3> t. The high 
temperature series therefore offers a quantitatively re- 
liable tool to compare to experiments. Our approach 
complements recent dynamical mean field studies that 
can be applied to lower temperatures 24( . We use up to 
10 th order in the expansion of the grand potential [2lj . 
£1 = — In (2)//3, to extract thermodynamic quantities for 
a uniform system (7 = 0). In terms of the series coeffi- 
cients, X( m \ the expansion reads: 



0Cl = lnz + ]T (Pt/z ) m xW(w, C), 
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where £ = exp(/3^) is the fugacity of the uniform system, 
Cl = fl/N, w = exp(-/3f7) and z = 1 + 2( + ( 2 w is 
the partition function of a single site in the atomic limit. 
In a local density approximation (LDA) we assume that 
each site of the trapped system can be approximated with 
parameters for a uniform system. With the replacement 
£ — > x t — cxp(/3^) the LDA becomes fi LDA = £\ CV 

To show that the LDA is an excellent approximation 
for the high temperature regime studied here we also 
compute the exact second order contribution to the grand 
potential in a trapped system: 
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where the second sum double counts nearest neighbors 
and 
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+ XiXj [I-S-U,8+U + h-U,-S+u] 



+ x 2 wls+u,-s-u + x 2 wI-s+uj~u- (4) 

In the above expression the quantities 5 = fii — fij 
and I A _ A = (exp(/3A) - 1 - /3A)/(/3A) 2 simplify to 
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Figure 1: Total compressibility (top panel), system radius 
(middle panel) and the core compressibility (bottom panel) 
plotted as a function of the chemical potential in the trap 
center for t — 0.054 kHz and U = 5 kHz. The arrows in the 
bottom panel indicate the Mott regime at T = 0.2 kHz. 



the uniform limit for 5 = and /q.o = 1/2. The 



LDA is recovered in the limit Xi 
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w)x 2 ./ (3U + x r (l + x 2 w) yielding a simple expression for 
the grand potential of a trapped system in the contin- 
uum approximation (valid for large particle numbers): 
-/?ft LDA » jd d r[lnz Q , r + z{t/3/z 0tr ) 2 X^ BA ], where z is 
the coordination number. This second order high tem- 
perature expansion holds for any bipartite lattice in any 
dimension. As we show below there are no discernible 
differences between the LDA and the exact second order 
results for the parameters studied here. 

We first focus on the compressibility which can distin- 
guish the incompressible ground state of a gapped phase 
from a compressible metallic phase in a homogeneous sys- 
tem. The compressibility per particle is defined as k = 
N-^idrii/d/i where N = - X> Ai (-/?0). 

We compute the compressibility with the 10 th order se- 
ries but find essentially no distinction from the second 
order results for (it < 0.9. The top panel of Fig. Q] plots 
the compressibility as a function of the chemical potential 
in the center of the trap, /1, for two different values of T. 
At low chemical potentials (—1.0 kHz < /i < 0.8 kHz) the 
system forms a compressible Fermi-lattice gas. At higher 
chemical potentials (0.8 kHz < /i < 4.2 kHz) the system 
forms a Mott insulator at the trap center. Here and in 
what follows we define the Mott regime with a central 
density (n*) = 1 ± 0.01. A weak, barely distinguishable 
feature around fi = 5.5 kHz in the top panel indicates 
that a finite fraction of the system has become incom- 
pressible. Otherwise, we find that n is a nearly smooth 
function of all experimentally relevant parameters: the 
total compressibility incorporates the edges of the sys- 
tem which overwhelms signatures of the Mott transition. 
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Figure 2: Compressibility (dotted line) and core compress- 
ibility (solid line) versus chemical potential in a uniform sys- 
tem, 7 = 0, with U = 7 kHz and t = 0.054 kHz. The top 
(bottom) panel sets T = 1 kHz (T = 0.5 kHz). The three 
incompressible regions correspond to a pinning of the density 
at 0, 1 and 2 for chemical potentials near fi — —U /2, U /2 and 
3(7/2, respectively. For fi < (7/2, n° vanishes because the 
system shows very little double occupancy. 



We next study the system radius, R 7 because the com- 
pressibility is not directly measurable in experiments. R 
is denned as the root mean square of the distance aver- 
aged with respect to the density [25[ . The central panel 
of Fig. [T] plots R in units of the Thomas- Fermi radius, 
Rtf = a(3A^/47r) 1 / 3 , where a is the lattice spacing. The 
plateau in the middle panel of Fig. [1] results from a com- 
bination of edge effects and incompressibility at the trap 
center. As the chemical potential increases, we add more 
particles to the edge gas, and the size R scales as the 
size Rtf, the relevant length scale for the edges. A simi- 
lar plateau can be seen in the compressibility in units of 
the Thomas- Fermi compressibility at a value k/ktf = 1- 
Nonetheless, we have again found only weak, edge de- 
pendent features that indicate the formation of a Mott 
insulator. 

We are thus looking for a robust technique that probes 
the gapped phase without resorting to edge effects. We 
will show that the core compressibility per particle fulfills 
this requirement: 



K C EE N^dD/djJL, 



(5) 



where the double occupancy is conventionally defined as 
D = 3(0,)/ dU = £,•<»<»<)■ Our definition of D dif- 
fers from the definition in Ref. [||. In Ref. Q a quan- 
tity related to its dimensionless derivative, 3D/dN, was 
measured, which directly relates to the ratio of core 
and total compressibility (without knowledge of /i) via 



3D/3N = (3D/dfx)(3n/dN) 



J /n. 



Figure 3: The core compressibility ratio versus on-site inter- 
action strength for \i = 4 kHz. The top panel fixes t = 0.054 
kHz for several temperatures, T — 0.07 kHz (dot-dashed), 
0.2 (dashed), 0.5 (solid) and 1.0 (dotted). The bottom panel 
fixes T = 0.2 kHz for t = 0.01 kHz (dashed) and 0.1 (dot- 
ted). The insets plot the density along a cross section in the 
cubic lattice for the same parameters as the bottom panel of 
Fig. [2] but in a trap and with /j, = 6.5 kHz (left) and fi = 3.5 
kHz (right). The circles are computed in a full second order 
expansion for the density while the solid line is computed in 
the second order LDA. The Mott gap pins the density near 
(m) — 1 at the trap center (right panel). For /i > (7/2 (left 
panel) a central compressible region arises with m > 1 at the 
trap center. 



In Fig. [2] we plot k c and k versus /i for two differ- 
ent temperatures in a uniform system and find that k 
is essentially identical to k when the system has doubly 
occupied sites, but is zero otherwise. This is a key fea- 
ture with useful implications for trapped systems: k c , 
by taking the derivative of the double occupancy, mea- 
sures the compressibility of the core region with density 
larger than one (see the insets of Fig. [3]) and is therefore 
insensitive to the edges. 

A comparison between the top and bottom panels of 
Fig. [2] shows that k and k agree at low temperatures 
and high chemical potentials. This can be understood by 
considering a single site in the atomic limit: (« c / ' n) t= o = 
(1 + C)[2 + C + W /C\^ 1 ■ At zero temperature we can see 
the implicit cutoff in the chemical potential (n c / n)t=o — > 
6(p — U/2). Measurements of 8D/dN therefore exclude 
low chemical potentials (and thus edge effects) in trapped 
systems, yielding a measure of the core compressibility 
ratio (CCR), k c / k, at low temperatures, T <C U. 

The LDA is an excellent approximation for the param- 
eters considered here. The circles in the insets of Fig. [3] 
plot the density computed without the LDA. Comparison 
with the LDA results (solid line) shows remarkable agree- 
ment for pt < 1. We have also compared the compress- 
ibility computed with and without the LDA and have 
found very little distinction (less than 4%) in the regime 



of validity of our high temperature series /3t < 1. The 
LDA is thus a good approximation and the high temper- 
ature series expansion yields a highly accurate tool for 
quantitative comparison with ongoing experiments with 
strongly interacting (t/U <C 1) fermions in optical lat- 
tices. 

We now make contact with recent experiments of 
Ref. @, reporting a measurement of the quantity Sd/5N 
where S stands for an approximate derivative taken by 
linear fitting to experimental data. This experiment 
measures the double occupancy fraction d = 2D/N in- 
stead of the double occupancy D. Estimating Sd/SN m 
dd/dN = {2k g /kN) - D/N 2 we find that the difference 
is small near unity filling implying that Ref. Q measures 
the CCR. Note that SN must be carefully chosen to en- 
sure the proper phase [3]. Fig. [4] compares 2k c / kN and 
Sd/SN, showing agreement near the Mott transition with 
only small differences of order D/N 2 . However, we as- 
sert that a measurement of the quantity SD/SN « k c /k 
provides a better and more direct connection to the core 
compressibility than Sd/SN . 

Both the CCR and Sd/SN plotted in Fig. H show a 
distinct signature of the transition to the incompressible 
Mott phase. At low U the core of the system lies in 
the compressible regime of the Fermi-Hubbard phase di- 
agram, i%i > 1. The peak structure in Fig. [4] originates 
from the choice of /j, and the peaks in Fig. [2] Upon in- 
creasing U we enter the Mott regime, = 1. (Note that 
compressibility alone cannot distinguish the Mott phase 
from a weakly compressible metallic phase.) Here the 
center of the trapped system opens a gap and the core 
compressibility drops exponentially to zero. The zeroing 
of the core compressibility (and the CCR) is therefore 
an indicator for the onset of a Mott insulating phase in 
the sample center. Recall that the total compressibility 
K shows very little structure as we enter the Mott phase 
of a trapped system, and the signal thus originates from 
the change in core compressibility (compare the top and 
bottom panels of Fig. [1]) . 

The inset of Fig. [4] makes a more direct comparison 
with Ref. Q . Here we match entropy in the dipole trap at 
a temperature Tjj with the entropy in the lattice at fixed 
N to obtain a lattice temperature T. The inset is plotted 
for the coldest temperatures ascertained in experiment. 
For larger temperatures, Tp /Tp — 0.195, where Tp is the 
Fermi temperature in the dipole trap, the Mott phase is 
less robust. 

We also find that the precipitous drop in core com- 
pressibility remains robust over a wide parameter range. 
The top panel of Fig. [3] plots n c / ' kN for several different 
temperatures. We find that the CCR shows a clear drop 
upon entering the Mott phase for temperatures well be- 
low U . For temperatures comparable to the Mott gap, 
the center of the trapped system becomes compressible 
and the clear signature of a Mott phase vanishes (dotted 
line). A temperature fixed point appears for U ~ fj, such 
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Figure 4: The core compressibility ratio (solid line) and the 
experimentally observed quantity Sd/SN (dashed line), plot- 
ted as a function of the on-site interaction strength for n = 4 
kHz, t = 0.054 kHz and T = 0.3 kHz. The central density is 
unity for U > 6 kHz. The two lines merge when there is little 
double occupancy. The inset plots the same but the values 
are obtained by matching the entropy of N = 80000 atoms in 
the dipole trap (at Td/Tf — 0.15) to those in the lattice with 
t = 0.054. 

that (H) « 0. Here one finds a crossover from a high 
temperature Fermi-gas to a Mott phase. The bottom 
panel of Fig. [3] varies the hopping to demonstrate that 
there is only a small shift with different hoppings. 

We have shown that observations of double occupancy 
of cold atoms in optical lattices reveal the core compress- 
ibility in a trapped Fermi-Hubbard model. This core 
compressibility clearly indicates the onset of incompress- 
ible states and describes recent measurements that show 
evidence for the Mott transition Q. Ongoing work will 
generalize our proposed technique to bosonic systems. 
The core compressibility implicitly excludes edge effects 
to reveal compressibility near the trap center. Transi- 
tions to incompressible phases (e.g., metal- insulator tran- 
sitions) nucleated at the sample center can be readily 
identified in experiments and compared with compress- 
ibility computed in uniform Hubbard models. 
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